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The structure function S(k; (3) for the one-dimensional one-component log-gas is the 
Fourier transform of the charge-charge, or equivalently the density-density, corre- 
lation function. We show that for \k\ < min (27rp, 27rp/3), S(k;f3) is simply related 
to an analytic function f(k; (5) and this function satisfies the functional equation 
f(k;P) = /(— 2/c//3; 4//3). It is conjectured that the coefficient of k 3 in the power 
series expansion of f{k\0) about k = is of the form of a polynomial in [3/2 of 
degree j divided by (/J/2)- 7 . The bulk of the paper is concerned with calculating 
these polynomials explicitly up to and including those of degree 9. It is remarked 
that the small k expansion of S(k; 0) for the two-dimensional one-component plasma 
shares some properties in common with those of the one-dimensional one-component 
log-gas, but these break down at order k s . 

Key words logarithmic potential; two-dimensional plasma; fractional statistics; random ma- 
trices; exact solution. 



1 Introduction 

The one-component log-gas, consisting of N unit charges on a circle of circumference length L 
interacting via the two-dimensional Coulomb potential $(r,r') = — log \ f — f'\, is specified by 
the Boltzmann factor 

A N,P II \e 2niXk/L - e 2 "^/Y, 0<Xj< L. (1.1) 
l<j<k<N 

1 Laboratoire associe au Centre National de la Recherche Scientifiquc - URA D0063 
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The constant Apf p, which plays no role in the calculation of distribution functions, results from 
scaling the radius of the circle out of the logarithmic potential, and also includes the particle- 
background and background-background interactions (a uniform neutralizing background is im- 
posed for thermodynamic stability). The thermodynamic limit N, L — > oo,N/L = p (fixed) is 
taken, which gives an infinite system on a straight line with particle density p. This system was 
first studied because of its relation to the theory of random matrices pCfl . The thermodynamic 
functions were obtained. The pressure P has the simple form 

(3P=[1- ((3/2)]p (1.2) 

at any inverse temperature (3. However, exact (simple) forms for the correlation functions were 
obtained by the pioneers only for the special temperatures corresponding to (3 = 1,2,4 (See 
Section 5). More recently, exact expressions for the two-body density were derived for arbitrary 
even integer Q and then for arbitrary rational f3 [11|. Unfortunately, these latter exact 



expressions are complicated multivariable integral representations which cannot be easily used 
as such for actual computations. The purpose of the present paper is to obtain explicit small k 
expansions for the structure function (the Fourier transform of the two-body density). 

The log-gas is an example of a system interacting via the (i-dimensional Coulomb system 
(here d = 2) but confined to a domain of dimension d—1. It therefore exhibits universal features 
— that is features independent of microscopic details such as any short range potential between 
charges or the number of charge species — characteristic of Coulomb systems in this setting 



[14]. One universal feature is the existence of an algebraic tail in the leading non-oscillatory 
term of the large-distance asymptotic expansion of the charge-charge correlation function. For 
general charged systems in their conductive phase, interacting via the two-dimensional Coulomb 
potential in a one-dimensional domain, this is predicted to have the form |jj 
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(1.3) 



/3(vrr) 2 ' 

where r is the distance. For the one-component log-gas, (11.31) can be verified for all (3 rational 



The verification is possible because the charge-charge correlation function (which for a one- 
component system is the same as the density-density correlation) is known explicitly for f3 
rational |ll[ (see Q2.3| ) below). In this work we further analyze the properties of the structure 
factor S(k; (3) (Fourier transform of the charge-charge correlation) for the one-component log-gas. 
In particular we are interested in the (3 dependence of the coefficients in the small k expansion 
of S(k;f3). 

The large distance behaviour ( |1.3| ) is equivalent to the small k behaviour 

S(* ;j S)~% (1-4) 

Furthermore, by making use of the equivalence of the charge-charge and density-density corre- 
lation in the one-component log-gas, together with the exact equation of state the second order 
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term in ( |1.4|) has been predicted for general (3 ||, giving 



Let 



f(k; (3) := ^-S(k; 0), 0<k< min (2vrp, vr/3p) 
| A; | 



(1.5) 



(1.6) 



and define / for k < by analytic continuation (we will see below that f(k; (3) is analytic 
for < \k\ < min (2-7T/9, 7T/3/9)). In Section 2 we use the exact result (|2.3| ) below to derive the 
functional equation 

The simplest structure consistent with ( |1.7| ) is 



(1.7) 



K(3 -S(k;P) = l + jr Pl (p/2)(B-y, \k\ < min(27rp,7r/3p) 



where Pj{x) is a polynomial of degree j which satisfies the functional relation 

Pj(i/x) = (-iyx~ j pj(x). 

Equivalently, ( |1.9| ) can be stated as requiring 

j 

Pj( x ) = J2 a iJ xl > a j,i = a j,j-i U even) 

1=0 

j-i 

pj(x) = (x-l)^2a jtl x l , ajj = a j:j ^i (jodd). 

1=0 



(1.9) 

(1.10) 
(1.11) 



Inspection of (1J3) shows that the conjectured structure (L8) is correct at order \k\ and 
furthermore gives 

Pl (x) = (x-l), (1.12) 



and thus a^o = 1 in (|l.ll|) . In Section 3 we use (|2.3|) to verify that the structure ( |l.q ) is correct 
at order k 2 and we compute P2(x) explicitly. In Section 4 we use an exact an exact evaluation 
of the two-particle distribution function for f3 even [H] to rederive the result of Section 3, and we 
also use this formula to verify the structure ( |1.8| ) at order k 4 and to compute p±(x) explicitly. 

Assuming the validity of ( |1.§| ) we see that Pj(x) can be computed from knowledge of the 
coefficient of \k\ 3 in S(k; (5), or the coefficient of \k\ 3 in d p S(k; (3)/d(3 p (p < j), for an appropriate 
number of distinct values of (3. Because the functional relation (|0]) has via ( |1.10| ) and ( |1.11 ) 
been made a feature of (L8) the values of (3 cannot be related by (3 t— > 4//3. In Section 5 the 
known exact evaluation of S(k; (3) to leading order in (3 is reviewed, as are the exact evaluations 
of S(k;2) and S(k;4). Also noted are the exact evaluations of S(k, 1) and S{k;(3) to leading 



order in 1//3, which according to (1.7) are related to S(k;A) and S{k;(3) to leading order in f3 
respectively. All of these exact evaluations are in terms of elementary functions, and so can be 
expanded to all orders in k. We then present the exact evaluation of dS(k; (3) j 'd(3 to leading 
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order in /?, as well as the exact evaluation of dS(k; (3)/dj3 evaluated at (3 = 2 and (3 = 4. The 
details of the latter two calculations are given in separate appendices. Again the final expressions 
can be expanded to high order in |fc|. Using this data all polynomials in the expansion ( |1.S| ) up 
to and including the term with j = 9 can be computed. This expansion is written out explicitly 
in the final section and some special features of the polynomials therein, relating to the sign of 
the coefficients and the zeros, are noted. A physical interpretation of the functional equation, 
based on an analogy with a quantum many body system, which identifies an equivalence between 
quasi-hole and quasi-particle states contributing to S(k;(3) for |A;| small enough is given. We 
end with some remarks on the possible occurence of a functional equation analogous to (|0]) in 
the two-dimensional one-component plasma. 



2 The functional equation 



The Boltzmann factor (1.1) also has the physical interpretation as the absolute value squared of 



the exact ground state wave function, |0) say, for the Calogero-Sutherland quantum many body 
Hamiltonian 

* = -£a£j + W2-i)(£) E sin27r(x . XhVL - 

j=l OX j L l<j<k<N Sm n \ X 3 Xk >l L 

This Hamiltonian describes quantum particles on a circle of circumference length L interacting 
via the inverse square of the distance between the particles. In the thermodynamic limit N, L — ► 
oo, N/L = p (fixed) the N particle system becomes an infinite system on a line with particle 
density p. The ground state dynamical density-density correlation function 

N 

/A n -(0, x; t) := (0\n(0)e- lHt n( X y Ht \0), n(y) := £ 5{y - Xj ) (2.2) 

3=1 



in the infinite system has been calculated exactly for all rational [11]. The fact that (|0)) is 



proportional to (LI) tells us that at t = ( p.2|) is equal to 

pj 2) (0,x) + p5(x), 

where pT^\ is the truncated two-body density, for the log-gas system. Thus the exact evaluation 
of 

S(k,P) := r (pT 2) (0,x)+ P 5(x))e ikx dx 

for the log-gas follows from the exact evaluation of (^^) for the quantum system. Taking /3 to 
be rational and setting 

P/2:=p/q=:X 

where p and q are relatively prime integers, the latter exact result gives || 

1 roo V 



1 poo V rl 

S(k; (3) = irC p , q (\) JJ / dr« II / d yj Ql q F(q,p, 
i=i Jo j=i Jo 



\\{ Xi , yj })6(k-Q Pjq ), (2.3) 
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where 



Cp,g(A) 

Qp,q 

F(q,p,X\{xi,yj}) 



r«(A)rp(i/A) 



2^ P \q\ nu r2 (p - a(» - 1)) n?=i r 2 (i - a - i)/a) 

9 P 

i=l j=l 

iw k» - ^'i 2A rw lyj ~ ^ v l 2/A 
nLinf =1 (^+A % ) 2 
i 

x n?=iteo* + ^)) 1_A ipS(i - %)) 1 - 1/A 



(2.4) 



In the domain of integration of (2.3) the integration variables are all positive and because of 
the delta function are restricted to the hyperplane 

q p 



\k\ 



3=1 



2vrp 



We see immediately from these constraints that the restriction yj < 1 in the domain of integration 
is redundant for 

\k\ < 2vrp. (2.5) 

Thus assuming ( [2.5| ) we can extend the integration over yj to the region (0, oo). Doing this and 
changing variables Xi i— > \k\xi and yj i— > \k\yj we see that for \k\ in the region ( |2.5| ) 

1 poo P roo 

S(k;P)=7T\k\C Ptq (X)l[ dxiH dy i Qj f ,F(g,p,A|{x i ,y i };A;)<y(l-Qp, ff ), (2.6) 

i=l J ° 1 = 1 ° 



where 



F(q,p, \\{xi,yj}; k) 



n?=ife(i + fcTiA)) 1 -* nf=i(%(i - fcy,-)) 1 - 1 ^ 

IW |a?i - ^i'| 2A IW \Vj - Vj'\ 2/X 



YiUxiU^ + xy^ 



Notice that (2/7) is such that the integral in ( |2.6j) is analytic for 

| A; | < min (2np, irpfi)- 
Thus according to the definition (|l.6| ) we read off that 



(2.7) 



(2.8) 



1 POO P POD 

f(k;f3) = 2ir 2 \C p , q (\)l[ d Xl \{ dy i Q; i ,F(g,p J A|{x <> %•};*;) 5(1 (2.9) 
i=i Jo i=i 

The functional equation (1.7) is a simple consequence of this exact formula. Thus we see 
that the integral in (2.9) is unchanged by the mapping A i— ► 1/A (and thus p «-» q) followed by 
k i— > —k/X. The precise functional equation (1.7) follows provided we can show that 

C m (X) = X 2 ^- 2 C q , p (l/X), 

which indeed readily follows from the definition of C Pi3 (A) in ( |2.4j ). 
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3 Expanding f(k\ (3) in terms of Dotsenko-Fateev type 
integrals 



Here we will develop a strategy based on the integral formula fl2.9| ) to expand f(k, (3) at order k 2 . 
This relies on our ability to compute certain generalizations of a limiting case of the Dotsenko- 
Fateev integral. This same method has been used in || |(| to compute the equivalent of f(k, (3) 
and its derivative at k = 0. 



We first expand the integrand in (|2.9| ) as a function of k. According to ( |2,7D we have 

oo 

F(q,p,X\{xi,yj};k) = G(q,p, \\{xi, yj}) (l + ^ H v (q,p, \\{x u y 3 })k v 



where 



G(q,p,X\{xi,yj}) 



rw \ x i - x v\ 2X rw \vj - yr\ 2/x 
nti w 3= m + \y 3 ? nti 4" A nu v) 



x+^hm.vA^W = n j.i(i + WA)'-»n5. 1 (i-*w)'-"»- {3A) 

The coefficients H v are homogeneous polynomials in {xi,yj} of degree za 
Let us now introduce the notation 

1 poo P poo 

Ip, q> x[h({xi,yj})} := JJ / d^iJJ / d?/,- Q 2 m G(q,p, X\{xi,y 3 }) 6(1 - Q P) q)h({xi,y 3 }). (3.2) 
i=i Jo i=l 

Because of the presence of the delta function the value of I p ,q : \ is unchanged if Q 2 q is replaced 
by Qp q for any n. Doing this and also introducing the usual integral representation of the delta 
function, we see by a change of variables as detailed in || that for h homogeneous of degree v 

T \ h (Sr ., 1M - J P,iXrM{xi,y 3 })} _ J P , q ,x[h({xi,y 3 })} 

h*M*m,Vji)\- (l/ + n _i)i - {y _ l )\ v- 6 ) 

where 

1 poo P poo 

Jp,q,\n[K{xhVj})] ■= II / dx iU / dyjQp.qGiQ'P'Mixi^j^e-^hdxuyj}) 
i=1 Jo j=1 Jo 

and Jp,q,X ■= Jp,q,X,0- 

Recalling (2.9) and (|3.1|) we see that in terms of the notation (|3.2j) 



f(k;P) = C Ptq (\)(lp iqt \[i\ + ^2lp t q t x[H u (q,p,X\{xi,yj})]k 1/ ^j. (3.4) 



The definition of H u in (3.1) shows 

H 2 ( q , P ,X\{^y 3 }) = (A 2 " A2 1)2 ( g« a - ^(|>? " A |>))> ( 3 - 5 ) 
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so to compute f(k,f3) at order k 2 our task is to evaluate 

q p 

I p ,g,x [Ql, q ] and I p ^x E A - A X] y f] ■ 

i=l i=l 



(3.6) 



Now because of the delta function in (3.2) 

^p,<j,A[Qp l(? ] = ^p,g,A [1] ? 

and we know from || that 



(3.7) 



C Pl ,(A)/ Pl9lA [l] = 1. (3.8) 

Thus our remaining task is to compute the second expression in ( |3.6| ) or equivalently, using (|3.3|) , 
to compute 

q p 

4 9 ,aE xj - X Vj] = lJp,qA X i} ~ X P J P,qAVj] ( 3 - 9 ) 

i=l j=l 

where the second equality, valid for any 1 < i < q and 1 < j < p, follows from the symmetry of 
the integrand. 

For this purpose we first note formulas for J p ^\[h] in the cases h = xf and h = y 2 . The 
formulas are 



2tt P 

2l (2g - 1/A + 1 



7rp 



+ Xyj 



Xq 



2vrp J P^,A[yi] np J P ,q,x 



Xi + Ay,- 



(3.10) 
(3.11) 



The derivation of ( 3.10| ) and (3.11) uses a technique based on the fundamental theorem of 
calculus. It was first used by Aomoto |IJ in the context of the Selberg integral, and has been 
adapted in (8| to the case of the Dotsenko-Fateev integral. 

Let us give the details of the derivation of fl3.10|) (the derivation of ( 3.1l| ) is similar). From 
the definition (13.11) we see that 



» c( , llh A|K»})-(^-2t^ + aA|:-i-)c(,, R *|{«,,»}). 



-A-l 



p[ Xi + Ay 



Thus 



" roo P roc Q . 

= IT / dxiJJ / dyj — (x 2 G(q,p,X\{xi,yj})e 
t=i Jo j=i Jo dXi 



-Q P 



P - x 2 



(A + 1) Jp, q ,\[ x i] ~ 2 X J P><?. A 



(A + l)Jp, q ,x[ x i] ~ 2 P J P,q,- 



Xi + Xyj 
2 



- X 2 



+ 2A E J p^ —^r - ^pJpaM] 



i'=l 



X 



Xi + Xyj 



+ 2X(q-l)J Pi 



Xi Xi' 

x 2 



Xi Xj/ 



2 ^pJp,q,\[x 2 i] 



(3.12) 
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where the first equality follows from the fundamental theorem of calculus, while the final equality, 
valid for any j = 1, . . . ,p and any i' = 1, . . . , q, {i' ^ i) follows by the symmetry of the integrand 
with respect to {x^} and {yj}- The symmetry of the integrand with respect to {x{\ also gives 



■I 



xj 



X 7 X ii f 



J, 



T 2 
.1 ;l 



X r) I X -) 



so we have 



J, 



Xi Xj> 



2 ( ^P>9> A 



x'i 



Xi Xj/ 



Xi' Xi 



Jp,q,\ ■ 



Substituting in ( 332 ) implies ( |3,10| ). 



From (3.10) and (3.11) we see that 
2i \„ 7 . r„.2i 

x P (2q - i/A + 1; 



q(2p - A + 1 



2tt P 

g(-A + l) 
2vrp 



<Jp,q,\ 



2tt P 



W-l/A + l) 



pq 

Tip 

A - 1 



Jp,q,\[yj] ^P,9,A 



xf - X 2 y 



(2np) 



Xi + Xyj 

2 Jp,q,\[Qp,q] = 



A- 1 



J p,qA l ] 
(3.13) 



Recalling Q, the results (|J), Q, Q and (|1|) give that 



I p ,q,\[H2(q,Pi M{xi},{yj}] 



1 (A - 1 



(2vrp) 



2 ^2 "W^M- 



Use of ( |3.8| ) then gives that the term proportional to A: 2 in ( |3.4[ ) is equal to 

It follows from this that the structure ( |1.8| ) is valid at order A: 2 with 

p 2 (aO = (x- l) 2 . 



(3.14) 



(3.15) 



4 Large-x expansion of pT 2 \(0,x) 

We have already remarked that the large- x expansion (|l.3| ) of the charge-charge correlation, or 
what is the same thing for the one-component log-gas, the large- x expansion of p^(0,x), is 
equivalent to the small-A: behaviour ( |1.4[ ) of S(k;0). More generally the expansion 

oo 
n=l 

is equivalent to the expansion 



s 



where the expansion ( f4.2[) contains the terms singular in k (i.e. of odd order in |fc|) only. This 
follows using the Fourier transform 



oo e ikx (—l) n \k\ 2n ~ 1 

dx = 7T- 

oo 



x 2n ■ f2n-l)! 



from the theory of generalized functions (see e.g. |]j 

From the equivalence between ( [4.1| ) and ( |4.2| ) we see the fact, following from ( 3.14] ) , that the 
term proportional to |/c| 3 in the small k expansion of S(k; 0) is equal to 

is equivalent to the statement that the term proportional to l/x 4 in the large x expansion of 
pj 2 ^(0,x) is equal to 

In this section we will derive (fD|) directly. We will also calculate the 0(l/x 6 ) term and so 
explicitly determine the 0(\k\ 5 ) term in (|4.2| ). 



The starting point for our calculation is an exact /^-dimensional integral formula for the 
two-particle distribution P(2)(®i x ) valid for even. With 

a ( , ^_ "tt 1 r(q + 1 + jc)T(b + 1 + jc)T(l + (j + l)c) 
n[ ' ' c) r(o + 6 + 2 + (AT +i - l)c)r(l + c) 

the formula gives that in the thermodynamic limit [Q] 

((/5/2)!) 3 e- 7ri/3 ^(27r / ox) /3 



P(2)(0,x) = p 2 (P/2f 



pi(3P/2)\Sf,(l -2/0,1 -2/0,2/0) 

( d Ul --- dup f[ e 2-P- Jn -l+2//3 (1 _ „ J-1+3//J -Q | Ufc _ |4//». (4.4) 



j=l j<k 



In a previous analysis [|| it has been shown that the non-oscillatory large-x behaviour is 
determined by the integrand in the vicinity of the endpoints and 1, with the requirement 
that 0/2 of the integration variables are in the vicinity of the endpoint 0, while the remaining 
0/2 integration variables are in the vicinity of the endpoint 1. Thus we write up/2+j = 1 — Vj 
(j = 1, . . . ,0/2) (this introduces a combinatorial factor choose 0/2 to account for the different 
ways of so partitioning the integration variables) and then expand the integrand (excluding the 
exponential factors which involve x) in terms of the "small" variables Uj, Vj (j = 1, . . . , 0/2). In 
particular we must expand 

13/2 /3/2 

nu - Uj r^(i - Vj )- i+ ^ n a - u i - ^) 4//3 . (4.5) 

3=1 l,l'=l 



The function (4.5) is a symmetric function of the variables {iij} and {vj} separately. Let {q K } K 
be a polynomial basis for symmetric functions with k denoting a partition (ordered set of non- 
negative integers) of no more than 0/2 parts, and suppose furthermore that q K is homogeneous 
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of order \k\ := k% + ■ ■ ■ + Kp/2- Then we can write 

(3/2 13/2 

U(i _ „.)-i+2aj (1 _ Vj) -i+w n (i - «i - 

3=1 l,l'=l 

= w ^Qk{ui, up^q^vx, v p / 2 ). (4.6) 

Substituting (|4.6j ) in ( [4.4|) , then following the procedure of ||, which involves extending 
the range of integration to Uj £ (0, oo), Vj G (0, oo) and changing variables Uj *— * 2iripxuj, 
Vj i— > —2mpxVj making use in the process of the fact that q K is homogeneous of degree |«|, we 
obtain the non-oscillatory terms in the large- £ asymptotic expansion of p^)(0,x). This reads 

P ^ X) ~ ^2)^S W - 2//?, I - 2//?, 2/0) 

re, v / 

where 

0/2 

i£/3 )/t := / diii • • • dup/2 ]J w^ 1+2//3 e _U! J| \u k - Uj^ 13 q K (u x , . . .,Up/ 2 ). (4.8) 

^[0,oo)5/2 J=1 j<fc 

The symmetry u> K)/i = evident from ( |4.6[) implies terms in (|4.7|) with |«| + \fj,\ odd cancel. 



Therefore the sum in (4.7) can be restricted to partitions such that |k| + \p\ is even, which means 
the asymptotic expansion only contains inverse even powers of x. 

To proceed further we must be able to compute the expansion coefficients w S)J1 as well as 
the integrals Kg r . For the former task it is convenient to choose q K equal to the monomial 
symmetric polynomial m K , which is defined as the symmetrization of the monomial x^ 1 • • • x^j^ 
normalized so that the coefficient of x* 1 ■ ■ ■ x^^ is unity. 

First, we have the well known expansion 

n 

Y[(l- Uj ) a = J2 a K m R (ui, . . . ,u n ) (4.9) 

3=1 l{n)<n 

where 

a K = n a k ■= -jf 1 ( 4 - 10 ) 

p=l 

with 1(k) denoting the length of k (i.e. number of non-zero parts). We can therefore immediately 
expand the first product in ( f4.6|) in terms of monomial symmetric polynomials. 

Consider next the expansion of the double product in (4.6). Making use of the formulas 

OO f \ 

(l- X f = ySZ^ljl x n (4.H) 

n oo 

II (£«**?) = E %- e ^a K m K ({ tj }) (4.12) 

3=1 k=0 £(K)<n 
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where a K is specified by the first equality in ( [4,10| ), we see that 



(3/2 

i[(l- Uj -v)W= £ (l-v) 2 -^c K m K ( Ul ,. 

3=1 £(K)<f3/2 



<U/3/2J 



(4.13) 



where 



P =i 



-4/i9)a 
fc! 



Expanding the factor (1 — v) 2 K l we can rewrite ( [4.13 ) as 

oo 

j - v flP = ^ w n (u!, . . . , u p/2 ;(3)v r 



(3/2 
3=1 



n=0 



for appropriate symmetric functions w n . Replacing v by vy and forming the product over j' 
using ( [4.12 ) we obtain 

(3/2 

II i 1 - U 3 ~ v j') 4/l3 = w a f w K m n{vi,...,Vi3/2) 

j,j'=l e(n)<(3/2 



where w K := ]lp=i w k v - The final step is to expand w$ /z l(K> w K in terms of {m^} and so obtain 
the expansion 



„0/2-*(k)„ 



(3/2 _^ 

II i l ~ u i _u i') 4//3 = ^2t^ K m ll (ui, ... ,Up/ 2 )m K {v u ... ,Vp/ 2 )- 
j,j'=i ma 



(4.14) 



The practical implementation of this procedure requires the use of computer algebra. We work 
with arbitrary (positive integer) values of (3/2. Furthermore, we only include terms with |/i| + 
\k\ < 6 throughout since according to (4.7) these terms suffice for the evaluation of the coefficients 
of l/x 2n , n < 3. 

Having obtained the coefficients t^ K in ( |4,14 ), we multiply the series ( PH ) with the two 
series of the form (4.9) representing the first two products in (4.5), expressing the answer in 
the form of ( |4.6| ), and so determining the coefficients w K ^. Again this step requires computer 
algebra. 



With w Ktfl in (4.6) determined, it remains to compute the multiple integral (|4.8| ) with 
= ma- For this task we introduce a further basis of symmetric functions, namely the Jack 



polynomials {P« (ui, • • • , Up/2)}- The Jack polynomials pj^^(zi, . . . , zn) with Zj := e 2mx ^ L , 
when muliplied by the ground state wave function |0), are the eigenfunctions of the Calogero- 



Sutherland Schrodinger operator (|2.1|) ||. Each polynomial is homogeneous of degree |k| and 
has the expansion 



pW( Zl ,...,z N ) 



m 



(4.15) 



where < is the dominance partial ordering for partitions: fi < k if |«| = \pt\ with k / /j and 
Ya=i fH — Sf=i K i f° r each p = 1, . . . , N. The coefficients a KM can be calculated by recurrence 



ii 



The significance of the Jack polynomial basis is that we have the explicit integral evaluation 
— II l dt l tte- t 'Pt ) (t 1 ,...,t N )Y[\t k -t j \ 2 / a = P^(l N )[a+(N-l)/a + l}t\ (4.16) 



W aaN 



j<k 



which is a limiting case of an integration formula due to Macdonald p9|| , Kadell [15] and Kaneko 
D). In ( p6| ) 



w aaN = n / d *i *? e ~*' n - ^ 

J=l*' 7<Jfc 



2/q 



AT-1 

n 

3=0 



r(l + (3 + l)/q)r(o + l+j/a) 
r(l + l/a) 



(a) 



N 

n 



~ (j ~ !)/« + K j) 
F(u - (j - l)/a) 



and Pi"'' (1^) denotes Pi"'' . . . , rcjv) evaluated at xi = • • • = re at = 1. 

To make use of ( f4.16| ) we must first express the monomial symmetric polynomials m K in 
terms of {P^ }^<k, which can be done using computer algebra from knowledge of the expansion 
( 4.151 ). Substituting in ( 4.161 ) allows the integrals K K p to be computed. 

After completing this procedure all terms in (4.7) for |k| + < 6 are known explicitly. 
Performing the sum and simplifying we obtain 



P(2)(0,x) ~ p I 1 



1 



+ 



309 



15(/3 - 2) 2 (/? 2 -3/3 + 4) 



+ 



(4.17) 



P(irpx) 2 2{3 3 (irpx)* 2(5 5 (tt P xY 
Note that this agrees with the known form ( jOj ) for the term 0(1 /x 2 ), and the form (f4.3|) for 
the term 0(l/x ). The term 0(l/x 6 ), due to the equivalence between ( [4.1| ) and (4.2), implies 
the term 0(|/c| 5 ) in the small-A: expansion of S(k,j3) is equal to 

3 



(/3/2-l) 2 ((/3/2) 2 -^(/3/2) + l)(i^) 5 



This is of the form of the conjecture ( |l.g| ) with 

Pa(x) = (x — l) 2 (x 2 



-x + 1). 



(4.18) 
(4.19) 



5 S(k] (3) for special /3 



Let us assume the validity of (1.8). The coefficients specifying the polynomials Pj(x) therein can 
be determined from knowledge of the coefficient of in S(k; (3) or d p S(k; f3)/d[3 p (p < j) at 

special values of (3. Now in the context of random matrix theory S(k; (3) has been evaluated in 
terms of elementary functions for /3 = 1, 2 and 4. The results are p0| 

' ?-gl°g(l + S), \k\<^P 

2-glog(^fe), |*|>2*P 
§, \k\<27rp 
1, > 2?rp 



5(A;;1) 



S(fc;2) 



(5.1) 



(5.2) 



5(fe;4) 



111 



M 

877 



log 



1 



2?rp 



|fc| < 47rp 
|fc| > 47rp 



(5.3) 
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Recalling the definition fll.6|) of f(k,/3) we read off 

mi) = i-^io g (i+— ) 



m-,2) 



1 — - log I I 



Tip' 

k 
2vrp 



(5.4) 
(5.5) 
(5.6) 



The exact evaluation ( |5.5|) implies that for all j Pj(x) contains a factor of (x — 1). In the 
case of j odd this gives no new information since the factor (x — 1) was already deduced as a 
consequence of the functional equation ( |1.9[ ). On the other hand, in the case j even this fact 
together with the functional equation ([O]) implies 



i-2 



Pj(x) = (x- lfj^bjjx 1 , 



(j even) . 



(5.7) 



1=0 



Consider now the constraints on the coefficients in (5.7) and ( |1 . 1 1| ) which follow from (5.4) 
and ( |5.6j ). As ( |5.4|) and ( |5.6| ) are related by the functional equation ( |1.7| ), and this is built into 
the structures (|5.7| ) and ( |l . 1 1 ) , only one of these exact evaluations gives distinct information on 
Pj(x). For definiteness consider (5.4). We see that 



(5.8) 



where the notation [k 3 \ denotes the coefficient of k 3 . Recalling (|1.§|), ( |5.7| ) and (1.11) this implies, 
for j even, 

j = 5 ((1 + 2-«- a ))6 i|0 + (2- 1 + 2-^ 3 %, 1 + • • • + (2"W + 2 -,72 )6 . + 

(5.9) 

while for j odd 
1 



J 



((1 + 2-^'- 1 ))a j , + (2- 1 + 2-^'- 2 ))a i , i + • • • + (2" 
+2-^- 1 )/^, _ 1)/2 ). 



■(j-l)/2+l 



+ 2-(- 1 )/ 2 - 1 )a J , (j _ 1)/2 _ 1 



(5.10) 



In the case j = 1 ( |5 . 1 Q|) gives S^q = 1 which reclaims ( 1.12j ), while in the case j = 2 (|5.9j ) gives 



•j,o 



1 which reclaims fl3.15|) . 



The exact form of S{k;(3) in the weak coupling scaling limit j3 — > 0, k — > 0, kj (3 fixed is also 
available. Introducing the dimensionless Fourier transforms 

/oo r oo 

(pL (0, x) + p6(xj) e ipxk dx, := p / $(x)e^' fc cix 

-oo ^ ' J —oo 

where <3?(x) := — log \x\ is the pair potential of the log-gas (thus the integral in the definition of 
the <3?(/c) is to be interpreted as a generalized function) we have [12] 

/9*(fc) 



S(k;0)~l 



(5.11) 
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Since 



7T 



and noting S(k; (3) = S(kp; (3)/ p we thus have that in the weak coupling scaling limit 



= p(l - — TjT 



+ |fc|/jr/3p/' 

Expanding ( |5.13| ) in the form ( |1.8[ ) and recalling ( [5.7D and ( 1.1 1] ) we deduce 

a,jfi = 1 and bjfi = 1 



(5.12) 



(5.13) 



(5.14) 



for all j. Using ( 5.14 ) in (|5.9| ) and ( 5.10 ) gives that in the case j = 3, cLp = — -g-, and in the 
case j = 4, bj t i = 
(|5.14|) gives 



f. The latter result reclaims ( |4.18| ) while the former result together with 



p 3 (x) = (x - 1)(1 - + x 2 ). 

O 



(5.15) 



An alternative way to derive (5.14) is to consider the (3 — > 00 low temperature limit. In 
this limit the system behaves like an harmonic crystal, for which we have available the analytic 
formula (|0 



p — — 00 



(5.16) 



where 



HP) 



1 /-V2 l - cos 2-irpt 
t-t 2 



vr 2 Jo 



dt. 



Taking the Fourier transform gives for \k\ < 2irp 

00 

S^(k;f3) = p J2 



e -k*f(p)/l3p 2 _ A e ikp/p 



p=—oo 

k 2 00 

' r p=— 00 



Jkp/p 



\k\M 
1 - |fc|/27rp' 



(5.17) 



This formula maps to the weak coupling result ( 5.1l| ) under the action of the functional equation 
(1.7) and so implies (|5,14|). 



6 Perturbation about (3 = 



The formula ( 5.13j ) is just the first term in a systematic weak coupling renormalized Mayer series 



expansion in (3. In the case of the two-dimensional one-component plasma, low order terms of 
this expansion have recently been analyzed by Kalinay et al. 0]. Results from that study can 
readily be transcribed to the case of the one-component log-gas. 



2 The denominator of the exponent in (3.10) of Q contains a spurious factor of it 2 which is corrected 

in 



(15-161 ) 
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Formally, the renormalized Mayer series expansion is for the dimensionless free energy (3F e * 
(in our (3F ex is written —f3F ex ), and one computes the direct correlation function via the 
functional differentiation formula 

c(0 X ) = ^ GX ) (6 l) 

The Ornstein-Zernicke relation gives that the dimensionless Fourier transform of the direct 
correlation function, c(k, (3) say, is related to the dimensionless structure function S(k; (3) by 

c(k;P) = l-^- — (6.2) 

V S(k;(3) ( 1 

so expanding c(k, (3) about (3 = with k/(3 fixed is equivalent to expanding S(k; (3) about = 
with k/(3 fixed. 



Now, transcribing the results of |L6| we read off that the weak coupling diagrammatic ex- 
pansion of c{x\,X2) starts as 

c(x!,x 2 ) = -P$(x 1 ,x 2 ) + ^(k( Xi ,X2)) 2 + ■ ■ ■ (6.3) 

where 



e ik(x 1 -x 2 ) 

K( Xl ,x 2 ) = -Pit I — — dk (6.4) 

3 \n\ ~r K 



with k = 0irp. This implies 



C{k](3) = ~]k\ + W-oo ^ \l\+K\ P k-l\+K (6 ' 5) 

The integral is straightforward (consider separately the ranges of I such that I > (I < 0) and 
pk — I > (pk — I < 0)). In terms of k' := pk/n = k/ir/3, 

c(t, (3) = -p + /? |fc/ | (2 + j fc | / j ;/|) log(l + \k'\) + 0(/3 2 ), (6.6) 
or equivalently using ([Tj) 

l + |fc/K| (1 + |A;/«|)(2 + \k/K\) 

Notice that the leading order term in ( |6.7j ) reproduces ( |5.11 ), 

The exact result ( |6.7l ) gives the explicit value of the coefficient of x in the polynomial Pj(x). 
Thus recalling ( l.ll[) and (5/7) we have 

l(b jA -2) = W\ {l + J )[2 + x) log(l + x), (J even) 

\{l-a hl ) = W\ {l + J ){2 + x) log(l + x), (J odd). (6.8) 

Furthermore, a simple calculation gives 
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so we have for example 
91 



05,1 



30 



31 

15 : 



07,1 



1607 
420 : 



263 
"~84 ! 



09,1 



791 
180' 



(6.10) 



Substituting 05,1 from ( |6. 10 ) and 05^ from ( 5.14 ) in (5.9) shows 05,2 = f§- Similarly, the 



value of &6,i above allows us to deduce that &6,2 = x - Thus we have 

91 o 62 91 



p 5 (x) 



(.-D(x 4 --x3 + -^-^ + l) 



37 



13 



37 



(x - 1) (x x 6 + —x 2 x + 1 



15 



15 



(6.11) 



We remark that according to the conjecture ( |1.8| ), the expansion of S(k,f3) about (3 = 
should have the structure 



S(fc, P) = /„(*/«) + + /^M*/*) + ■ • • 



(6.12) 



where 



= u j (cj t o + c^iu H ). (6.13) 

Consideration of the analysis of [16| reveals that the structure ( 6.12j ) will indeed result from the 
weak coupling expansion, however the structure ( |6.13 ) is not immediately evident. (Of course 
the explicit form f 2 as revealed by (|6.7[) exhibits this structure.) 



7 Perturbation about (3 = 2 and (3 = 4 

A feature of the couplings (3 = 1, 2 and 4 is that the n-particle distribution functions are known 
for each n = 2, 3, ... ]20|j . Introducing the dimensionless distribution 

£f(£l, ...,£„):= /)(„)(xi, . . .,X n )/p n 

we can use our knowledge of g(xi, . . . , x n ) for n = 2, 3 and 4 at these specific /3 to expand 
g(xi, X2) about /3 = fta to first order in /? — /3q. Thus with <3?(xi, 22) := — log \xi — x 2 \ we have 
II 



g(xi,x 2 ;P) = g{xi,x 2 ) + {/3 - /?o)j - g{xi,x 2 )$(xi,x 2 ) 

-2p / (5(^1, x 2 , x 3 ) -g{x 1 ,x 2 )) $(xi,x 3 ) dx 3 

J— 00 

1 /"°° / 

~t;P 2 / (^(^l) ^2, s 3 , x 4 ) x 2 )5(x3,X4) -g(xi,x 2 , x 3 ) 
£ J —00 ^ 

-#(xi, £2,^4) + 25t(xi,x 2 ))$(x3,x 4 ) dx 3 (ix4| + 0((/? - /3 ) 2 ) (7.1) 

where on the right hand side the dimensionless distributions are evaluated at (3 = (3q. Here we 
will compute this first order correction, and the corresponding first order correction for S(k; (3), 
in the cases (3q = 2 and (3q = 4 (we do not consider /? = 1 because of its relation to /?o = 4 via 
the functional equation (|Q|)). 
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j,k=l,...,n 



Now, in the case (3q = 2 we have 

g{xi,...,x n ) = det P 2 {xj,x k ) 
while in the case (5q = 4 

. . . ,x n ) = qdet P^Xj^Xk) 

where 

P A (Xj,X k ) = 



sin7rp(x — y) 
vrp(x - y) 



j,k=l,...,n 



sin 2-Kpxjk 



2npx jk 
1 <i ^sin27T / ox :; fc 



Si(27T/9Xjfe) 

sin 2irpXjk 



(7.2) 



(7.3) 



(7.4) 



27r/? dxjfc V 2-irpXjk 1 2irpxj k 
with Xj k := Xj — x k and Si(x) denoting the complimentary sine integral, defined in terms of the 
sine integral si(x) by 

00 sini 

at = — + si [x), si[x) := — I 

o 



r% sin t 7r /* c 

Si (x) = / dt = — |- si (x), si (x) := — 

Jo t 2 J x 



I x t 

In (^) qdet denotes quaternion determinant, which can be defined as 

l 

= V (-l) n - r 

j,k=\,...,n 



dt. 



(7.5) 



qdet P^Xj^Xk) 



E (- i r~' II ( P ^a, X b )Pi(x b , X C ) ■ ■ ■ P 4 (x d , x a j) 

Pes n i 



(0) 



(7.6) 



where the superscript (0) denotes the operation |Tr, P is any permutation of the indicies 
(1, . . . ,n) consisting of / exclusive cycles of the form (o b c- ■ ■ — > d — > a) and (— V) n is 
equal to the parity of P. Note that this reproduces the definition of an ordinary determinant in 
the case that P4 is a multiple of the identity. 



The task now is to substitute (7.2) in the case /3q = 2 and (7.3) in the case (3q = 4, and to 
compute the integrals. Consider first the case (3q = 2. After some calculation (see Appendix A) 
we find 

1 /sin7rpx\2 sin27r/?x 



g(0,x;P) 



1 



+ ci (2irpx) 



2(irpx) 

where C denotes Euler's constant while 



irpx ' l 'Z \ irpx ' 2irpx 

+ 1 2 ((log 2tt P \x\ + C) cos 2irpx - ci (2vrpx)) | + 0((J3 - 2) 2 ) (7.7) 



■ , \ ^ , , rcosi-l , 

ci(x) = G + log |x| + / dt 

Jo t 



cos t 



dt 



(7.8) 



denotes the cosine integral. From this we can compute (again see Appendix A) that up to terms 
0((/3-2) 2 ) 



jg + (/3 _2)pUlog(l 



S(k;(3) 



(2tt P ) 



p+( / S-2)pUlog 



1 1^^. \k\+2np \k 
\k\-2irp ~ 4irp 



log( 



1 



(27TP) 2 
k 2 



|fc| > 27T/9. 



(7.9) 



Let us consider the consequence of ( |7,9|) in regards to the expansion fll.q ). For \k\ < 2irp 
we observe that all terms but the one proportional to \k\ are even in k. This is consistent with 
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Pj(x) having the quadratic factor (x — l) 2 for j odd (recall ( |5,7|) ), but only a linear factor for j 
even (recall Ql.llp ). Moreover, we can use (7.9) to derive a linear equation for the coefficients 
{cij}. First we differentiate ( |7.9| ) with respect to /3, set (3 = 2 and expand about k = to obtain 



dS(k;f3) 



dp 



(3=2 



2 2vrp Z - 



1 / |*| \2jf 

j= J 2j(2j - 1) V^J ' 



|*| < 27rp. 



Recalling (|0|) and ( |5.7[ ) this in turn implies 

1 l/ _ „ \ 

2j(2j - 1) ~~ 2 V i -1 ' + 2a2 i- 1 . 1 h 202J-1J-2 + a 2 j-ij_ij. 

In the case j = 4 we deduce from this equation, ( 5.14|) , (|6.1[l| ) and ( |5.10|) that 



, , , 1607 2011 2 911 s 2011 4 1607 = B 
p 7 (x) = (x - 1)( 1 + — x - tt^x + -^rx t^ttX* + x° 



420 



280 



420 



(7.10) 



(7.11) 



Consider now the case (5q = 4. Due to P\ in (7.3) being a 2 x 2 matrix, the calculation 
required to compute (|7.1| ) is more lengthy and tedius than in the case (3q = 2, although the 
the common structure of n-point distributions means the two cases are analogous. Some details 
are given in Appendix B. Our final expression for g(x±, X2', /3) is given by flB.4p , We find its 
Fourier transform can be computed explicitly in terms of elementary functions, together with 
the dilogarithm 



dilog (x) :- 



logt 



dt. 



h l-t 

Explicitly, with p = 1 for notational convenience, up to terms 0(((3 — 4) 2 ) 



(7.12) 



7T 



S(k, (3) = S(k, 4) + ((3 - 4) ( - — + B (k) + 2B x {k) - 4B 3 (k) + 2B 5 (k) + B 6 (k) - B 7 (k)) (7.13) 



where 
B (k) 



3 3|*| \k\ , /4vr + |*k /_ 1 . _ 2 

■2 + ir + 47 log (Tr) + ( c+ 2 log(16 " 

I A; 
16 

2vr - \k\ 



\k\ 
2^ 



- (dilog ( 
log 



1*1 



2vr + |* 
_\k[ 
2vr 



dilog ( 



47T + |*| 
27T+ |*| 

*| < 47T, 



log 



1 



2tt 



log ( 



47T + |*| 



(7.14) 



S (*) 



-dilog 



|*| +4vr 
|*| -4vr 
|*| 



+ 



|*| , /* 2 - 16vr" 

|*| + 47T 



JM. 

16vr 



1 dilog ( • 



|*| 



* + 2vr 



dilog (t 



Bi(k) 



|*| -2irJ °V|*| + 2vr 

7T A IM i |fe| 1„„ 



dilog ( 



i _ 1*1 

1 2tt 



0, 



1*1 - 47T 

|*| - 2vr 

|*| < 47T 
1*1 > 47T, 



|*| > 4vr,(7.15) 
(7.16) 



B 3 (k) 



3 3 * 

h - 1 — 1 + 

2 8vr 



c(i 



|*| 
2^ 



1 3|*K, 
1 1 8 " 32^) bg 



2vr 
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B 3 (k) 



B 6 (k) 



+^(log|l - H|) 2 + ^(4vr - |*|) log(4vr - |*|) - -L|fc| log \k\ + ilo g 4vr 



+ 



64vr 
\k\ 



(dilog ( 



2vr 
1*1 



B 6 (k) = 



32vr V 
+21og(27r)log 
0, \k\ > 4vr, 

B 5 (fc) 



_3 3jfc| 
2 8vr 
J*|_/r[ 
32vr V 3 



+ 



7T 



2vr+ \k\J ' 12 
|*| 



i - 



2vr 



+ log(27r + |A;|)log|l-^|+5 2 (fc)), |*| < 4vr, 




M 

2?r 



|*| > 47T 



16vr 



log 



1 - 



\k\ 
2vr 



+ 



(C + l0g(47T-|fc|))(l - 



\k\ 1*1 

'47T - \k\ 



1 - 



(7.17) 
(7.18) 



\k[ 
2vr 



-2 log 



|*| 
2vr 

0, |*|>4vr, 
S 7 (fc) = 



log( 



4tt - |*| 
2^ 



+ g±{k) , |*|<4tt, 



o, 



IAI 

2?r 



) 2 , |*|<4vr 

|*| > 47T, 



(7.19) 
(7.20) 



2?r 



2tt- 



£ - log (l-Jg) log 1^1 < 2vr 



with 
51 (fc) = 

52 (fc) = 
S3(fc) = 

54 (A;) = 

{ dilog(^j-2dilog^^j+log^^jlog^^_j, 2vr < < 4vr. 

(7.21) 

The above formula for S(k; 0) in the case \k\ < 2ir (recall here p = 1) can be used to expand 
dS(k; f3)/d(5 about k = 0. For this task we use computer algebra, which gives the result 

5fc 4 3|*| 5 27fc 6 



dilo g(£qf) " dilog (, 

dilog(^_) - dilog( w ^_) - £ + log (^) log (^L) , 2vr < |*| < 4vr, 
' dilog(^g)-^+log(2vr-|A:|)log(l-g), |*| < 2n 
^ -dilog(^_) + log(4vr - |*|) log (M - l) , 2tt < |*| < 4n 
~ 2dilog(ft) - log (g) log \k\ < 2n, 

K dil °g(^) - d H%) + lo g(^) lo g(S))' 2vr< |*| <4vr, 



' -*-dilog(^)+2dilog(^)+21og(^)log(^), |*| < 2. 
, dilo g(^)- 2dUo g(wF^) +1 °g(^) lo g(A)' 2vr < < 4vr. 



dS{k;(3) 



dp 



(3=4, 



\k\ |fc| 3 

16vr 256vr 3 3072vr 4 



+ 



37|*| 7 



+ 



1273* 8 



245760vr 7 18350080tt 8 
1949|*| n 



+ 



4096vr 5 81920vr 6 

887|*| 9 4423* 10 



+ 



27525120vr 9 293601280vr 10 



2752512007T 11 



+ ••• 



(7.22) 
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This allows us to deduce a further equation for {bs,j}j=o,...,4 and {a9j}j=o,...,4> which in combi- 
nation with fl7l0D , ( |6~T0| ), ( gig ), ( pT9l) and ( p0|) implies 



p 8 (x) 



(x 
(x 



l) 2 (l-^ 

7 V 84 
791 



1697 



11 



x + 
84 315 

73603 



180 X+ 7500 



6337 

1008 
7355 

~504 



x° + 



■x° + 



1697 

315~ 
2231 

135' 



263 



-x 



84 
7355 

"504 



x° + x" 



(7.23) 



■x J + 



73603 
7560 ^ 



791 
180 



x' + x l 



(7.24) 



8 Conclusion 



Collecting together the evaluations fli~T2D , plo] ), fTT5D , ( |4T9D , ( pH"D , ( f7TTTD , ( [7T23D and ( [^24]) , 
and substituting in ( |1.8j ) we have that for \k\ < min (2np, Trfip) 

7T/3 



1*1 



■S(fc;0) = 
1 

+(x - l)y 
+(x-l)V 

+(x — l)(x 2 X 

6 



11 +i)y 3 



+(x-i)V-§ 

91 3 



- \x + 1)/ 



I : " - 62 n 91 ^. r 

; 4 x d H x 2 x + l)y 5 

15 30 



' lH - r 30-' .., 
+(x-l) 2 (x 4 -|> + ^x 2 -|x + l) y 6 



' v 15 4 15 

w fi 1607 = 2011 4 911 , 2011 2 1607 , 7 

+(x - l)(x b x 5 H x 4 x 3 H x 2 x + l)y 7 

v A 420 280 105 280 420 ;tf 

x2/ fi 263 „ 1697 4 6337 , 1697 , 263 . „ 

+(x - 1) (x 6 x 5 H x 4 x 3 H x 2 x + l)u 8 

J y 84 315 1008 315 84 

1W o 791 7 73603 fi 7355 = 2231 4 

- l)(x 8 x 7 H x 6 x 5 H x 4 - 

/v ISO 75R0 504 1 35 



+0(y 



,10% 



84 

7355 , 73603 , 791 „, „ 

x d H x z x + l)y y 

504 7560 180 1 

(8.1) 



where x = (3/2 and y = \k\/iT(3p. With the coefficient of y 3 denoted Pj(x) as has been throughout, 
we recall from our workings above that po(x), pi(x), P2(x) and Pa{x) have been calculated for 
general values of (3. In all other cases the calculation has relied on the assumption that the Pj(x) 
are indeed polynomials. On this point we remark that in such cases, excluding j = 8 and 9, we 
have more data points than is necessary to uniquely specify Pj(x), assuming it is a polynomial, 



and our extra data points are consistent with the explicit forms presented in (8.1). 
We remark that the structure exhibited by 



l|) is familiar from the study of exactly solv- 



able two-dimensional lattice models [10|. In this field one encounters two- variable generating 
functions G(x, y) say with series expansions of the form 



G(x,y) 



n=0 



x)y ' 



1.2) 
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in which H n (x) is a rational function, and furthermore the denominator polynomial in H n {x) 
only has a small number of (typically no more than two) distinct zeros. For example, the 
two-dimensional Ising model with couplings J\ (J 2) between bonds in the horizontal (vertical) 
direction and x := exp(— AJi/ksT), y := exp(— 4J2/&;bT) has for its spontateous magnetisation 
the celebrated exact expression (see e.g. ||) 

When written in the form (18.21) one finds 



where P n (x) is a polynomial of degree 2n — 2 which satisfies the functional relation 

P n (x) = x 2n - 2 P n (l/x). (8.5) 



As emphasized in [10], the exact solution fl3.3|) can be uniquely determined by the functional 
form ( |8.4[) , together with the functional (inversion) relation ( |8.5[) and the symmetry relation 
M(x,y) = M(y,x). For the structure function of the log-gas we have no analogue of the 
symmetry relation and so cannot characterize ( |2.6[ ) this way. 

One immediate feature of the polynomials Pj(x) in ( |S.1| ) is that for j even the polynomial 
Pj(—x) has all coefficients positive, while for j odd the polynomial pj(—x) has all coefficients 
negative. Another general feature of the Pj(x) in ( |8.1| ), obtained from numerical computation, is 
that all the zeros lie on the unit circle in the complex x-plane. This can be rigorously determined 
numerically because the symmetry ( |1.9[) implies that if xq is a zero of Pj(x), then so is 1/xq, 
which will be the complex conjugate of xq if and only if IxqI = 1. 



The quantum many body interpretation of (1.1) allows us to give a physical interpretation to 



the functional relation (|1J) . As the functional relation is derived from the integral representation 



), it is appropriate to recall [11] the physical interpretation of that formula. In ( |2.9[) , with 
(3/2 = p/q, there are q integrals over x\ € (0, 00) and p integrals over yj 6 (0, 1). The variables 
Xi can be interpreted as being rapidities of quasi-particle excitations, while the yj are rapidities 
of quasi-hole excitations. Thus the transformation (3 1— > 4/(3 is equivalent to interchanging p and 
q and thus the quasi-holes and quasi-particles. In (^) this does not lead to an integral of the 
same functional form as before; although the functional form of the integrand is conserved, apart 
from a renormalization of k, the domain of integration is different for {x^} and {yj}- But with 
k restricted as in (^) both sets of variables can take any value in (0, 00). The quasi-particles 
and quasi-holes play an identical role and the functional equation results. 

It is of interest to consider the small k expansion of S(k; T), T := q 2 jk^T (q = charge), for 
the two-dimensional one-component plasma. As mentioned earlier, this has recently been the 



object of study of Kalinay et al. |16 |. They obtain results which imply 



2 " r 5(Jfe;r) = 1 + (- - 1)— + (- - -)(- - -){—) 2 + 0(k 6 ). (8.6) 



k 2 y ' ' v 4 '2vrrp v 4 2 /v 4 3'\2irTp 
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where k := |fc|. The structure of ( |8,6|) bears a striking resemblence to flS,l| ) with T/4 correspond- 
ing to x and k 2 /2i:Tp to y. In particular with g(x,y) := (27rr//c 2 )S'(A:; T), the expansion (8.6) 
to the given order is such that 

g(x,y) = g(-;-yx). (8.7) 

Furthermore, writing 

oo 

g(x,y) = l + J2Mx)y l (8-8) 
i=i 

we have u\(x) = (x — 1), U2(x) = (x — 3/2) (x — 2/3) so uj(x) is a monic jth degree polynomial 
for j < 2. However we can demonstrate that this analogy breaks down for the I = 3 term in 



To demonstrate this fact, suppose instead that the functional equation ( |8.8[) was valid at 
order / = 3 in (|8.Sj) and «3(x) is a monic polynomial. Then 113 must be of the form 



u 3 (x) = (x-l)(x 2 + ax + l). (8.9) 
From the definition of g(x,y) we can check that this is equivalent to the statement that 

I /7TT>x4 r r 8 s(r . r)cF=(4!) 2 (x _ 1)(x 2 + ax + -g_ (810) 
P V 2 / JR2 



But as noted in [16], it follows from the perturbation expansion of [13| that 



K^) 4 L r8s(r;r)df - -* + (r-2)4!(£^-2) + o((r-2) 



k=0 

-4! + (r-2)4!^ + 0((r-2) 2 ). (8.11) 



The term in (8.11) proportional to T — 2 is incompatible with ( |8.10| ) which gives instead 



1 8 

(T - 2)4! T 



independent of the value of a. Indeed in |1^] evidence is presented which indicates us(x) is an 
infinite series in x, although we have no way of determining if the functional equation (^^) also 
breaks down at this order. 
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Appendix A 

In this appendix some details of the derivation of ( |7.7| ) and ( |7.9| ) will be given. To simplify 
notation we take p = 1 throughout. The first step is to substitute ( fT% ) and (|7.1|) and simplify 
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by expanding out the determinant and cancelling terms where possible. This shows that up to 
terms O((0 - 2) 2 ) 



= 1 - (P 2 {xi,x 2 )f + 03-2){ - (1 - (P 2 (x u x 2 )) 2 )$(x u x 2 ) 

(■Z 

-2 



(- (P2(X 2 ,X 3 )) 2 ~ (P 2 (X!,X 3 )) 2 + 2P 2 (x 1 ,X 2 )P 2 (x 2 ,X 3 )P 2 (x 3 ,X 1 ))$(x 1 ,X 3 ) dx 3 

1 f°° ( 

- \ \^P 2 {x\,x 3 )P 2 (x 3 ,x^)P 2 (x^x\) - 4:P 2 (xi,x 2 )P 2 (x 2 , x 3 )P 2 (x 3 , x A )P 2 (x Al xi) 

^ J — oo ^ 



-2P 2 (xi, x 3 )P 2 (x 3 , x 2 )P 2 (x 2 , x 4 )P 2 (x 4: , x{) + 2\^P 2 (xi,x 3 )j yP 2 (x 2 ,x 4: )j )$(x 3 , x 4 ) dx 3 dx A j. 

(A.l) 

The convolution structure 

/■oo 

f(yi - x)g(x - y 2 ) dx 



often occurs in the above integrals. Such an integral can be transformed by introducing the 
Fourier transforms / (<?) according to the formula 

l r°° 

f( Vl - x)g(x -y 2 )dx = — f(l)g(l)e- il(m - y2) dl. (A.2) 
Making use of this formula typically leads to simplifications. 



For example, consider the first integral in ( A.l ). Starting with the Fourier transform 

00 sin 2 ™ ( 1 - jg, \k\ < 2tt 



-e tkx dx = { zn (A.3) 
W | 0, \k\ > 2vr 

and ( 5.12| ), application of ( |A.2| ) gives 

/oo , . 2 / , 2i" , \k\\ 7T dk 

[P 2 (x 2 , x 3 )j ^(x 3 , Xl )dx 3 = J ^ (l - — cos/cxis — . (A.4) 

This expression is indeed simpler than the original, but it suffers from being ill-defined, due 
to the singularity at the origin. However its derivative is well-defined, and can furthermore be 
evaluated in terms of elementary functions giving 

d . . . sin27rx 1 . . 

Tx A ^ x) = ^r-x (A ' 5) 



Also, we have 



00 sin 2 ttx 
ix 

oo 



A l (0) = - dx- — r=-log|sI = C + log27r-l, (A.6) 

{7TX) Z 



where C denotes Euler's constant. Together ( |A.5[) and ( |A.6| ) imply 



At (x) = - S ™ 2lTX + ci (2ttx) - log \x\ (A. 7) 

2ttx 



where ci(aj) denotes the cosine integral (|7j). 
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The other six integrals in ( A.l ) yield to similar techniques. We find 
A 2 := I (P 2 {x 1 ,x 3 )) 2 <f>{x 1 ,x 3 )dx 3 = A x (0) 

P 2 (x 2 ,x 3 )P 2 (x 3 ,xi)^(x 3 ,x 1 )dx 3 = J [C + log(7r + k)J coskx 12 — 

1 / . .\sin7rx12 1/./ \ vr\cos7rx 
= - C + log2vr - log 212 + ci(27txi 2 ) - si(27rxi 2 ) + - 

2 V / 7rxi2 2 v 2 1 vrxi2 

/oo 
P 2 (xi,x 3 )P 2 (x3,X4)P 2 (x4,xi)$(x3,X4) dx 3 dx^ = A 3 (0) = Ai(0), 
-oo 

/oo 
^2(^2, x 3 )P 2 (x 3 ,X4)P 2 (x4, ^1)^(^3, 0:4)^x3^x4 = A 3 (x 12 ), 
-00 

/oo 
P2(xi, x 3 )P 2 (x 3 , x 2 )P 2 (x 2 , X4)P 2 (x 4 , xi)$(x 3 , x 4 ) dx 3 dx A 
-00 



2(irx 12 ) 



— -2 (C + log 2ir + cos 27rxi2 ( log |xi 2 1 - ci (27rxi 2 )^ - sin 27rxi 2 (si (27rxi 2 ) + — 



(P 2 (X1,X 3 )J (P 2 (x2,X 4 )J $(x 3 ,X 4 )(iX3dx4 = J ^ — (l - —j — COS fcX!2 

1-COS27TX12 sin27TXi 2 
= -log X12 ^ hci(2vrxi2), (A 

(2TTX 12 y 27TX12 

where si(x) denotes the sine integral defined in (|7 



Of the results (|A.8|) , the evaluation of Aq is the most difficult, so it is appropriate to give 
details in that case also. We observe that Aq consists of the convolution of P 2 (xi, X3)P2(x3, X2) 
regarded as a function of X3, and <I>(x3,x 4 ), and P2(x 4 , xi)P 2 (x2, x 4 ) regarded as a function 
of x 4 . It simplifies the calculation to take as the origin in both integrations the centre of 
the interval between particle 1 and particle 2, which is achieved by the change of variables 
X3 !— > X3 + (xi + x 2 )/2, x 4 1— > x 4 + (xi + x 2 )/2. Use of (|A.2j) then shows 



/oo ^t, . \ 2 TT 

^sin^7T- ¥f)x 12 , \k\ < 2tt 

oo 7TX13 ^X 32 Q, > 2vr 

where the second equality in ( |A.10|) follows after further use of ( |A.2| ). Thus 

J-2tt 27T V 7TX / \k\ 

As in ( |A.ll ), this integrand is ill-defined. To proceed further, we write 

,4 6 (x) = 4 1) (x)+4 2) (x) 



\rn \ f°° sin7rxi 3 sin7rx32 , , I nx 12 Dil1 V" 2 , -. . . 

X12) := / cosKX3ax3 = < (A. 10) 



where 



(i) f 00 dA; f /sin(7r — \k\/2)x\ 2 /sin7rx\2 , l ir 

d 2tt\\ vrx / \ 7tx / Jjfcj 



(2) _ /sinvrxy dfc vr 

^ ^ " V^^J J-2n2^]k\' 
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The integral defining A^ is well defined and can be computed by elementary means. The 
integral defining is singular. It coincides with the singular part of -Ai(O) (recall (A. 4)), and 
so from (A.6) we have 



Sin 7TX\ 2 



TTX 



(C7 + log27r) 



Collecting together the above evaluations of A1-A7 and substituting as appropriate in (|A.1|) 
gives ( ff7|) . 

The next task is to evaluate the Fourier transform. Now the evaluations of A\ and A7 are 
given as Fourier integrals, so their Fourier transform is immediate: 

7T 1 



FT Ai (x) 
FT^ 7 (x) 



Ifcl 2 



, \k\ < 2ir 



7T Ifcl 



(A.12) 



while for \k\ > 2ir 

FTAi(x) = FTA 7 (i) = 0. (A.13) 

We can check that the constants A2 and A4 cancel when substituted in (A.l), and so play no 
further part in the calculation. 

Of the remaining terms, consider first the first term proportional to — 2 in ( |A.l| ), Aq(x) 
say. Making use of (|A.2j) we see that 



7T 



FTA (x) = -—+ — 



2tt 



dl TC 



\l\ 

\k\ ' 7-2tt 2vr \l - k\ V ~ 2vr 
For \k\ < 2tt minor manipulation allows the singular part 

r2w dl TV 



2tt 2TT 



C + log 2vr 



(A.14) 



(A.15) 



to be separated, while the remaining convergent integrals are elementary. We thus find that for 
\k\ < 2vr 



FT M*) = ~ + {C + log 2. + I log (l - (A) 2 ) } (1 



1*1 
2^ 



WW ^ /*+]*! 

2tt 2tt s V fc 



)■ 



(A.16) 



For > 2ir the integrals in ( |A.14| ) are convergent and also elementary. In this case we find 



1 



4vr z 



™ . , s 1 , Ifcl + 2vr Ifcl 



log ( 



(A.17) 



To compute the Fourier transform of Aq, we begin by making use of (|A.2| ) in ( A. 9 ) thereby 
obtaining 

FT,4 6 = 



where 



V(l, k) := / dx V(l, x)e ikx = X\k\<-n-\l\/ 



25 



with the equality in the latter formula following from the explicit form ( A,10| ) of V(l,x) and 
then computation of the resulting integral, and where \T = 1 for T true and xt = otherwise. 
Thus 



oo 



, dl f 00 dki tt 

FTA 6 (X) = I ^ — J ^ — X|fexl<vr-|i|/2X[fcx-fe|<vr-|Z|/2 jjj 

(l _ M) (C + log 2vr) - (l - S) log ^ - ^(2vr - |*|), |*| < 2vr 
0, |*| > 2vr 

(A.18) 

where use has been made of the generalized integral evaluation ( [A. 15 ). 

The final Fourier transform to consider is 

^^sinvrx . , . ^^sinvrx f n dk\ 1, , , . 1, , , ,\ 

FT A 3 (x) = FT / _i(C + -log(7r + *i + - log tt - *i ))e tkix 

ttx irx J-tt 2ir V 2 2 / 

1 f°° r 1 1 

= 2vr 7 ^ Xle[-n,7T]Xle[-7T+k,7T+k] [C + - log(vr + I - k) + - log(vr - (Z - *)) 



where to obtain the equality use has been made of (|A.2j ). Evaluating the integral gives 

FT ^ 3 <*> = (c +1 o g2 ,) (l - - H log fi + i(, - H) log(l - M) (A . 19) 

for |/c| < 27r, while for [*| > 2ir 

„ m sin-/rx . . 

FT A 3 (x)=0. A.20) 

7TX 



Substituting the above results as appropriate in the Fourier transform of ( |A.1| ) gives the 
result (171). 



Appendix B 



In this appendix we outline some details of the calculation of (7.1) in the case /3q = 4 and 



show how this leads to ( [7.13 ). Because J7.2|) and (|7.3| ) formally have the same structure upon 



expansion (recall the definition of qdet (|7.6[) ), the formula ( |A.l| ) formally maintains its structure 
when generalized to the case /3q = 4. Thus we have 

g 2 (xi,x 2 ;f3) 

= 1 - (P 4 (X1, X 2 )P 4 (x 2 , X{j) + (P ~ 4) | - (l - (P 4 (xi, X 2 )P 4 (x 2 , Xl)) (0) )$(xi, X 2 ) 

/oo . 
( - (P 4 (x 2 ,X 3 )P 4 (x 3 ,X 2 ))(°) - (P 4 (X 1 ,X 3 )P 4 (X 3 ,X 1 ))^ 
-oo 

+2 (P 4 (a?i , x 2 )P 4 (x 2 , £3 )P 4 (x 3 , si ) ) (°) ) $ {x x , x 3 ) (fx 3 

'Y POO , 

~o / (4(^4(3:1, x 3 )P 4 (x 3 , x 4 )P 4 (x 4 , xi)) (0) -4(P 4 (xi,x 2 )P 4 (x 2 ,X3)P4(a;3,a; 4 )P4(x 4 ,xi)) (0) 
2 7— 00 ^ 

-2(P 4 (xi, x 3 )P 4 (x 3 , x 2 )P 4 (x 2 , x 4 )P 4 (x 4 , Xi)) (0) 

+2(P 4 (xi, x 3 )P 4 (x 3 , xi)) (0) (P 4 (x 2 , x 4 )P 4 (x 4 , x 2 )) (0) ) *(x 3 , a*) dx 3 dx 4 ] + 0((/3 - 4) 2 ). (B.l) 
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We treat each of the seven distinct integrals in (|B,1| ) in an analogous way to their counterparts 
in (A..1), although extra working is involved due to P4 being a matrix rather than a scalar. 

The final results are 



Bi{x 12 ) 



B> 

B 3 (x 12 ) 



B 5 (x 12 ) 



B 6 (xi 2 ) 



B 7 (x 12 ) 



(P±{x 2 , a;3)-P4(^3, ^2)) ( $(X1,X3) dx 





-00 

^ dk 

-4tt 27T 



\k\ \k\. 



cos kxi 2 



log|xi 2 | 



\k\ 

sin47rxi2 , cos 2irx 12 

Cl (47TX12J Dl(z7TXl2j 



47TX12 



47TX12 



— 00 

DC 



(0) 



(Pi(xi, a; 3 )P4(x3, xi)J ' '3>(xi, x 3 ) (fx 3 
Pi(x 2 , x 3 )Pi(x 3 , xi)^(xi, x 3 ) dx 3 



Bi(0) 



C + log4vr- - 



" |/ 3 (X1 2 ) -|/o Xl2 (/lW + /2(*))* 

-3/1(^12) + 3/3(^12) 3/1(^12) + \hix 12 ) 
f 00 ( \(o) 

/ [P4{xi,X3)P4 : {x3,xa)Pa{X4,xi)j $(xi , x 3 ) dx 3 dx A = Bi(0) 
^4(^2, 2:3)^4(^3, x 4 )P 4 (x 4 , xi)$(x 3 , x 4 ) dx 3 dx 4 



00 
00 



|/l(xi 2 ) + |/ 2 (X1 2 ) " |/ 3 (X1 2 ) -imfl(t) + \f 2 {t) ~ \f 3 {t))dt 

-(3/1(^12) + 1/2(^12) - i/^xia)) |/i(x 12 ) + f/ 2 (x 12 ) - |/3(^i2) 

(0). 



[Piixi, x 3 )Pi(x 3 , x 2 )-P 4 (x 2 , x 4 )P 4 (x 4 , x 3 )J $(x 3 , x 4 ) dx 3 dx A 



4 2pfe| ((^(k'^ 12 )) 2 + cos(A;xi2/2)5i(A;,xi2)(— g 2 {k,x u ) - —g 3 (k,x 12 ) 



+ cos(kxi 2 ) i—g 2 {k,xi 2 ) 
+ ^^ 9 s(k,x 12 ))( 



ik 
8^ 



33(fc,«12) 



47T 4-7T 

sin(|fc|xi 2 /2) 



4tt 



52(^,^12) 



(4vr 



-cos((2tt- |fc|/2))xi 2 ) 



4tt 



Attx\ 2 



sin((2vr - |fe|/2))x 12 ) 



2vrx 2 2 



r 00 / \ (0) / \ (0) 

J [Pa{xi,x 3 )P4(x 3 ,xi)) [PA{x 2 - l x A )P A (x A ,x 2 )j $(x 3 ,x 4 ) dx 3 dx 4 



In 



1 



— + — In 

Air 8ir 



2vr 



+ 



log [x| 

Itt 
- Itt 



sin 47rx 



47TX 



+ ci (4-7rx) 



7r e?fc 
R COS ^ 12 2^ 
Si (2-7tx) sin 27rx Si(2-7rx) cos 2-7rx 



8tt 2 x 2 



Attx 



\k\ log 



\k\ 2 , dA; 

COS kx rr 

2vr 128vr 2 



(B.2) 



where 



fx(x) := / (2C7 + ln(4vr 2 -k 2 ))coskx^- 

J-2n 27T 

sm27rx,„ , ,,. cos27rx r ,. /- . 
-(C + ln47r - In |x| + ci(47rx)) Si(4ttx) 



7TX 



7TX 



27 



h{x) 

gi(k,x) 
g 2 (k,x) 
gz(k,x) 



-— Si(2vrx), f 3 (x) :-- 

7TX 

sin((27r - |fcl/2)a?) 
2ttx 

ci((2vr - |Jfc|)as) - ci (2nx) 



sin 2ttx 



TTX 



Si((27r 



)x) + Si (2ttx) 



(B.3) 



When substituted in ( |B.1| ), the constant terms B 2 and B\ cancel, and we obtain the formula 

/ \(o) 
g(xi,x 2 ;P) = 1 - \Pa{xi,x 2 )Pa{x2,x{}\ 

+ (/? " 4)|( - 1 + (P 4 (X 1 ,X2)P4(X2,X 1 )Y ^(X 1 ,X 2 ) 

+2B 1 {x 12 ) - ^P 4 (x 12 )B 3 (x 12 )) {0) + 2(p 4 (xi 2 ) j B 5 (x 12 )) ( ° ) 
+B 6 (x 12 ) - £7^12) j + 0((/3 - 4) 2 ). (B.4) 



We will demonstrate the close analogy with the (3 = 2 calculation of Appendix A by giving 
the derivation of the integral formula in (|B.2j ) for Bx(x± 2 ). As in the derivation of the integral 
formula ( A.4j) for Ai(x\ 2 ), our strategy is to use the convolution formula (AJ2). However here 
the Fourier transform of Pa(xi, x 2 )P l i(x 2 , x{) is not immediate. What is immediate is the Fourier 



transform of Pa(xi,x 2 ). Thus from the definition (7.4) we see that 

/oo 
Pi{xi,x 2 )e Lkxu dx l2 
-00 



1/2 i/2k 

-ik/2 1/2 





\k\ < 2vr 



\k\ > 2ir. 



Use of (|A.2|) then shows that for \k\ < Att 
FT P i (x 1 ,x 2 )P 4: (x 2 ,x 1 ) := 



2vr 
-2tt 



1/2 i/2l 
-il/2 1/2 



1 \M i Ifcl 1 



1 _ JM 



1/2 -i/2(k-l) 
i(k-l)/2 1/2 





dl 

X\k-i\<2n 2~ 



l_M + Mi c 

- 1 4vr T 8tt iu e 



1 _ JM 



while for Ifcl > 47T 



(B.5) 



(B.6) 



FTP 4 (^i,x 2 )P 4 (x 2 ,xi) = 0. 

The results ( |B.5[) and ( |B.6[) are the analogue of ( |A.3 ) in the working leading to the evaluation 
of A\{xi 2 ). The integral formula for B\{x\ 2 ) in ( [B.2D now follows from ( |B.5| ), ( p3.6| ) and ( p. 12 ) 



upon a further application of ( A. 2 ). 

The Fourier transform of flB.4j ) can be computed explicitly. The final result has already been 
stated in ( [7.13 ). This is obtained through the intermediate results 



FTBj(x) = Bj(k) for j = 0, 1,3,5,6,7 



28 



with 



(o) 



-B0O12) := [Pa(xi,X2)Pa(x2,xi)J $(x 1 ,x 2 ) 

and the Bj specified by ( |7,14| )-( |7.20| ). We will illustrate the working by giving some details of 
the computation of Bo(k) for \k\ < Air. 

Using and (|U2|) we see from (]Aj) that 



FT Bo (^12) = /_ 47r + 11 log 

(" 



1 - li. 
2vr 



7T dZ 



~ 2^ 



7T dZ 



fc - Z| 2vr 

47r 7T dZ_ 

-4tt I A; — Z| 2tt 



\k\ 
'8^ 





2-7T — |Z| 


/ log 


2tt- 


J -in 



fc-/| 2tt 

7T dl 



2tt 



7T dZ 



\k - l\ 2ir 



I 2ir 



(B.7) 



where the second equality, which follows from minor manipulation of the first integral, is moti- 



vated by the desire to separate the singular integral. Thus in the second equality of (B.7) only 
the first integral is singular. It is essentially the same as the first singular integral in ( A.14| ), 
and is evaluated as 

r 4 ^ 7T dl 



(B.8) 



The second integral in the second equality of flB.7| ) also appears in the evaluation of ( A.14| ). 
An elementary calculation shows 
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To evaluate the third integral in ( |B.7| ) we suppose without loss of generality that k > and 
write 
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dl 



The only non-elementary integral is the second term of the first integral. This can be computed 
by checking from the definition ( 7.12| ) that for —An < I < 



4 [ dilog f — -) + log 1 + — log f — 
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^k + 2-K 
In total we therefore have 
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To evaluate the final integral in ( B.7 ), a similar approach to that leading to the evaluation 
(B.12) is adopted. Minor complications arise because of the need to modify the formula (B.ll) 
for I > k. We find 
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where gi is defined in (|7.21|). Substituting ([eO|)-([B.13D in gives the result ([7.14[). 
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